Generators of nonclassical states by combination of the linear coupling of boson 
modes, Kerr nonlinearity and the strong linear losses 
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We show that the generators of quantum states of hght can be built by employing the Kerr 
nonlinearity, a strong linear absorption or losses and the linear coupling of optical modes. Our 
setup can be realized, for instance, with the use of the optical fiber technology. We consider in 
detail the simplest cases of three and four coupled modes, where a strongly lossy mode is linearly 
coupled to other linear and nonlinear modes. In the three-mode design, our scheme emulates the 
third-order nonlinear absorption, allowing for generation of the single photon states, or the two- 
photon absorption allowing to generate the phase states. In the four-mode design, the scheme 
emulates a non-local absorption which produces an entangled state of two uncoupled modes. We 
also note that in the latter case and in the case of the phase states generation the output state is in 
the linear modes, which prevents its subsequent degradation by strong losses accompanying a strong 
Kerr nonlinearity. 
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I. INTRODUCTION 



Robust, efficient and reliable generators of nonclassi- 
cal states of light are needed for the current and possi- 
ble future technology. Nonclassical states are the work- 
ing tools of quantum information and computation [1], 
besides other applications in quantum diagnostics and 
tomography, metrology, communications, etc. An ideal 
quantum state generator can be thought as a device pro- 
ducing a given quantum state by pushing a button, i.e. 
deterministically. It is known that a nonlinear media 
can be used to produce a nonclassical state of light from 
a classical (i.e. coherent) one [2], for instance, the n- 
photon Fock state can be obtained, at least in principle, 
in this way [3]. However, the experimental realization of 
the deterministic nonclassical state generators faces such 
great difficulties that the "determinism" requirement is 
currently dropped (it is pertinent to note that the overall 
efficiency of one of the most promising candidates for the 
"on-demand" single photons generation, the color centers 
in diamond, does not exceed a few percent [4]). Instead, 
one has to be currently satisfied with the conditionally 
generated nonclassical states. For instance, the genera- 
tion of the single photon states [2, 5] is currently based 
on producing the twin photon pairs by the parametric 
down-conversion or the heralded single photons gener- 
ated in the mesoscopic pn junctions [6, 7]. 

The main difficulty of the deterministic nonclassical 
states generation using a nonlinear media is the losses, 
since a strong nonlinearity required for such a genera- 
tor is inevitably accompanied by a strong absorption (in 
particular, this is the case of the nonlinear glasses used 
for fabrication of the nonlinear optical fibers). On the 
other hand, tailoring of a classical state to a nonclassical 
one by using the nonlinear absorption was also consid- 
ered. Thus the dissipation, considered usually as a de- 
structive effect to the fragile quantum features (e. g. the 
so-called "Schrodinger cat" quantum state is a prover- 
bial example of this [8]), can be used also to drive the 



system into a nonclassical state (this is the essence of the 
"quantum state protection" [9]). A wide variety of initial 
mixed states (or even an arbitrary mixed state) can be 
transformed into a desired pure state in this way, thus 
providing an extreme robustness of the scheme [10]. The 
artificially designed dissipative systems were also shown 
to be the means of the quantum computation [11]. 

The problem is that the practical realization of a given 
nonlinear loss channel for photons is rather difficult. One 
can engineer the dissipation needed for generating of the 
nonclassical states more or less freely working with the 
ions in a magnetic trap [12, 13], or with the atoms in the 
optical lattices [14]. However, designing similar reservoirs 
for photons is a much more challenging task, where even 
the theoretical schemes are mostly lacking. A notable 
exception is a couple of schemes proposed recently for 
generating of the phase states (JO) + e''^\l))/V2 [15-17] 
and the single-photon states [18]. Alas, they still require 
the very specific nonlinear dissipation (for example, the 
third-order nonlinear absorption is needed for the scheme 
of Ref. [18]). Remarkably, in a recent work [19] we have 
shown that one can emulate the needed nonlinear absorp- 
tion by a simple design based on the optical coupler. 

The purpose of the present work is to study in a more 
general setting the potential of our method of the quan- 
tum state generation, employed in Ref. [19]. The central 
idea behind our scheme is the emulation of the needed 
nonlinear absorption by using a combination of the lin- 
ear absorption, the Kerr-type nonlinearity and the mode 
coupling. It can be explained as follows. Consider a 
strongly lossy local mode linearly coupled to several other 
local modes. On the short evolution time (or propagation 
length), inversely proportional to the strong absorption, 
the lossy mode is emptied, whereas the rest of the modes 
are still significantly populated (curiously, this is due to 
the quantum Zeno effect [20]). Then the evolution can 
again be represented as that of a single lossy mode cou- 
pled to other modes but now in a coherent basis (i.e. a 
superposition of the local modes). Since some of the local 
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modes are nonlinear - the usual Kerr nonlincarity is as- 
sumed - the coiipling of this new lossy mode to other co- 
herent modes is now nonlinear. When the coherent lossy 
mode is also emptied, the dynamics of the remaining co- 
herent modes is governed by the nonlinear absorption, 
which can be tailored by selecting the particular setup of 
mode coupling. We show that, additionally to the single- 
photon states, also the phase states and the entangled 
states can be generated by our method. Moreover, out- 
put states in the second and third aforementioned exam- 
ples are localized in the linear modes of the system, thus 
preventing the degradation of the output quantum state 
by the strong losses, always accompanying a strong Kerr 
nonlinearity. 

The paper is organized as follows. In section II we 
describe the general scheme for the quantum state gen- 
eration. We focus on the three coupled modes in section 
III, where we consider two nonlinear side modes, subsec- 
tion III A, and one nonlinear and one linear side modes, 
subsection III B, coupled to a common lossy mode. More- 
over, we study the generation of multi-photon Fock states 
by this scheme in subsection III C. In section IV the case 
of four coupled modes is considered. In the Conclusion 
we outline the perspectives of the method. 



II. THE GENERAL SCHEME FOR 
EMULATION OF NONLINEAR LOSSES 

Our design for emulation of nonlinear absorption in- 
volves a boson mode subject to strong linear losses and 
coupled to other linear and nonlinear boson modes. We 
assume that the dynamics can be described by the mas- 
ter equation in the standard Lindblad form [21]. For the 
system with Hamiltonian H and a strongly lossy mode 
h, with the loss rate F, the evolution of an arbitrary op- 
erator A{t) is given by the equation 



h: 



HA 

^=i[H,A]+TV[b^]A, 



(1) 



where t is the evolution variable which will be called 
"time" for below (we use dimensionless variables), 
V[b^]A = b^Ab - ^{b^b,A} (here and below the capi- 
tal letters denote the time-dependent operators, while 
the small letters are used for the time-independent bo- 
son mode operators). We assume that the Hamiltonian 
H can be expanded in the powers of the lossy mode as 
follows 



H= J2 H<^°''^\b'')"b^. 



(2) 



a, 13=0 



Here the coefficients H^"'^^ are operators in all other 

modes of the system. We limit our consideration to the 
Kerr nonlinearity, thus a,/3 <2 (note, however, that our 
method works for arbitrary nonlinearity). 

Any time-dependent observable A{t) can always be ex- 
panded in the infinite power series with respect to b^ and 



A{t) = J2 A^'"'^Ht){b^rb'^ , 

a,0>O 



(3) 



where the coefficients A^°''^\t) arc time-dependent func- 
tions of the creation and annihilation operators of the 
other modes of the system. 

For strong losses, i.e. F ':$> for the evolution 

times t 3> 1/F the mode 6 is almost empty and one can 
adiabatically eliminate it from the operator equation (1). 
The key observation is that the leading order term in the 
expansion, i. e. A^^'^\t), is coupled by the Hamiltonian 
of Eq. (2) to just a few higher-order terms. The latter 
can be easily found in the explicit form in the adiabatic 
approximation. Inserting Eqs. (2) and (3) into Eq. (1) 
we have 



dAm 
dt 



^[^(0.0)^ ^(0,0)] _ ,^(0,1)^(1.0) ^ ,^(0,1)^(1 



,0) 



+2ii?(0'2)^(2,0)_ 2^^(0.2)^(2.0)^ 



(4) 



Observing that for large times, t^T, the term A^"''^) (t), 
with a -|- ^ > 0, decays as exp{— F(a -I- (3)t/2} we re- 
tain in the evolution equation for this term the only non- 
decaying term and the term multiplied by a large param- 
eter (the decay rate F) 



; W I 

dt 

which results in 



[if(",/3),A(0,0)]_£(^_t^^(a,/3)^ (5) 



[if("''^),A(0'0)(i)], a+/3>0. (6) 



F(a + /3) 



Finally, inserting the result (6) into Eq. (4), we arrive 
at the master equation with generally nonlinear artificial 
losses 

Q=1.2 

where the Lindblad generator is defined below 

Eq. (1). Eq. (7) has the decay rate in the usual form 

of the quantum Zeno effect, i.e. the actual losses are 
inversely proportional to the bare loss rate F [20]. 

Let us make the following observations about Eq. (7). 
Though we have assumed that the terms and 
fj{0,2) jji ii^Q expansion (2) of the Hamiltonian involve 
the nonlinearities, the strongly lossy mode is not required 
to be nonlinear (in fact, due to the strong loss, its non- 
linearity is always negligible). Moreover, it is only lin- 
early coupled to other modes of the system. Since Kerr 
nonlinearities are always diagonal in the spatially local 
boson basis, one needs the coupling to more than one 
mode, which makes such coupling diagonal in the non- 
local coherent basis and allows to emulate some nonlinear 
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absorption. Thus, a hncar couphng of a strongly lossy 
mode to two or more other modes, at least one of which 
is nonlinear, leads in the second adiabatic reduction to 
emulation of a nonlinear (coherent) lossy mode b coupled 
to other nonlinear (coherent) modes. 

We should mention the conditions which must be im- 
posed for the adiabatic elimination of the lossy modes 
(in the two stages of the scheme) and the more important 
condition of the scheme efficiency related to the existence 
of the natural or induced (see below) linear losses. These 
conditions are particular to each realization of the general 
scheme and arc given below for each example. We only 
mention that the first type of conditions are not actu- 
ally necessary but only sufficient for the emulation of the 
nonlinear absorptions, but they allow an analytical treat- 
ment via the adiabatic elimination, as outlined above. 
The second type is more important, since, as we will see 
below, the linear losses, even not naturally present in the 
output modes, appear due to coupling to the lossy modes 
though with a significantly reduced rate (as in the case of 
the phase state generation, section IIIB, and the genera- 
tion of the entangled state, section IV). Moreover, both 
types of conditions point to the conclusion that a strongly 
asymmetric modal coupling is needed for the design of a 
specific nonlinear dissipation and the localization of the 
output nonclassical state in the local modes. 

Finally, the system described by Eq. (7) can be real- 
ized using a fiber coupler or multi-core nonlinear fiber, 
where the optical mode of the central core is linearly 
coupled to the optical modes of the side cores/fibers, 
where some of the latter are made of the nonlinear media 
[19]. In this case, the cores have to be the single-mode 
ones, however, the consideration can be generalized to 
the multi-mode case as well. Another possibility is to 
implement the electromagnetically induced transparency 
(EIT) media, where one can achieve a very high Kerr 
nonlincarity [22] (termed even the "giant Kerr nonlin- 
earity" ) . Co-directional modes travelling through such a 
media can interact simultaneously and resonantly with 
the strongly damped emitters, in this way realizing the 
correlated loss of Eq. (2). A variation of this scheme is 
the propagation of several differently polarized modes in 
a nonlinear fiber with a strongly dissipative impurities 
(atoms, quantum dots, etc.) present in the core. When 
the average distance between impurities is much larger 
then the pulse length, the correlated modal loss can be 
achieved. 



III. NONCLASSICAL STATE GENERATION BY 
THE THREE MODE COUPLER 



cases and the possibility of localization of the generated 
state in the mode with low losses. We consider the fol- 
lowing Hamiltonian 



{gidi +g2a2)al + h.c 



2) 

(8) 

where the strong loss is i n mode oq with the loss rate Fg 
such that Fq G = v's'i + ^^'^ ^ I '^1,2 1 (here 
are real). The master equation for an observable 
A = A{t, aj, ao, a\,ai, a2, a2) of the whole system reads 



^=i[H„A]+J2^,V[a]]A. 



(9) 



The first reduction of Eq. (9), given by Eq. (7) with 
the expansion in Eq. (3) in powers of aj and ao, features 
ij(o.i) = Y, 9ja\ and il^^'^) = 0, and emulates the lossy 

i=l,2 

coherent mode hi (with the linear losses rate 4G^/ro) 
coupled to another coherent mode 62- The new modes are 
given by the following rotation of local modes {q = (71/52) 




(10) 



Eq. (7) now involves the Lindblad terms for modes 61,2 
and the Hamiltonian describing the interaction of modes 
hi and 62, which are derived from Eq. (9): 



j=l,2 



E 



g^F, 



+ 



g(r2 



h]Ahr--{h]hy,A} 



(11) 



^'°'°^= E Tf-2(9&]• + (-l)^■^'l')'(9&i + (-l)^■^'.•')^ 

(12) 

where j' = 2 for j = 1 and vice versa. 

We assume that the parameters are chosen in such way 
that mode hi has the shortest evolution time and can be 
adiabatically eliminated. We consider two particular lim- 
iting cases. In the first one both modes 01^2 are carried 
by the nonlinear media with the same nonlinearity, thus 
Ki = K2. In the second case only one of the two modes 
is nonlinear. As we shall see, there are significant differ- 
ences between these cases. 



Let us see how the above ideas can be applied for 
the case of three interacting modes. Recently, we have 
considered some particular realizations of this system 
[19, 23]. Here we analyze in detail the possibilities of the 
three-mode arrangement for generation of the nonclassi- 
cal states, discussing the symmetrical and asymmetrical 



A. Two nonlinear side modes: generation of the 
single-photon state 

Let us consider first the case of two nonlinear side 
modes 01,2 , when one has ki = K2 = k and Fi = F2 = F. 
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FIG. 1. (Color online) (a) Dynamics of poo (dashed line), pii 
(solid line), p22 (dash-dotted line) and average photon number 
(dotted line)as given by the solution of Eq. (13); F = O.I73. 
(b) Dynamics of pn as given by Eq. (13) for the linear loss 
rate F = O.OI73, O.I73, 73 (solid, dashed and dash-dotted 
lines). Dotted lines depict the solutions given by Eq. (19) 
for the corresponding values of F. (c) Dynamics of the ratio 
P22/P11 for F = O.I73, 73, IO73 (thick solid, thin dashed, 
thin solid lines). Thick dashed, dash-dotted and dotted lines 
correspond to the same values of the linear rate, but with 73 
being a hundred times smaller, (d) Dynamics of P22 as given 
by Eq. (13) for the linear loss rate F = O.I73 and F = 73 
(solid and dashed). Dotted lines depict the solutions given 
by Eq. (19) for the corresponding values of F. In all figures 
72 = O.OI73, and the solution to Eq. (13) is found for the 
initial coherent state with the amplitude z = 7. 



This case was considered before [19] as the scheme for 
the single-photon state generation. Here we turn our at- 
tention to the important problem of interplay between 
the hnear and nonlinear losses arising on the long evolu- 
tion time scale. Also, we outline some general features of 
dynamics produced by the artificial nonlinear losses. 

In Ref. [19] it was shown that after the modes ao and 
bi are emptied, the leading-order master equation for an 
observable A — A{t, b\,h2) becomes 



^ «i[i/(o-"),A]+r2?[6^]A- 
with the decay rates 



(13) 



72 



(1 



r)' 



73 



(1 



and the Hamiltonian 



_H-(o,o) 



1) 



2 I "2. 



(14) 



(15) 



2(i+g2)2 ^ v(i+'?': 

We note that the symmetric coupling, i.e. q = 1 leads 
to 73 = (see Eq. (14)) and was considered recently in 



Ref. [23] in the context of cold bosons in the triple-well 
trap. One of the notable results is that the two-photon 
dissipation (with the rate 72) cannot be captured by the 
mean-field approach, usually applied to systems of large 
number of bosons, though the induced linear dissipation 
is described by the mean-field model. Thus, quite nat- 
urally, one can expect of such a dissipation to generate 
nonclassical states. 

When the linear and two-photon losses are absent 
(r = and 72 = 0), the scheme described by Eq. (13) is 
able to generate the single-photon state from a large am- 
plitude coherent input [18, 19]. Indeed, the purity of the 
output state can be easily estimated. The three-photon 
dissipation (with the rate 73) leaves the following opera- 
tors invariant: ^00 = |0)(0|, ^01 = |0)(1| and ^10 = A\^. 
Since for large times the density matrix p involves only 
the Fock states |0) and |1), the following relations are 
valid: 



Poo(oo) = poo(O), pio(oo) = pio(O), pii((X)) = ^pfefe(O). 

fc=i 

(16) 

Eq. (16) gives for the coherent input state p(0) = \z){z\ 



the following estimates: poo ~ e '^1' 



, Pw 



"1^1 and 



Pii 



Thus, for |z| > 1 and i > I/73 the den- 
sity matrix p is exponentially close to the single photon 
state |1)(1|, i.e. 1 - Tr{p2} « 2e-l^l'. 

Since neither of gj can be zero, the two-photon losses 
are also inevitable. However, for a strong asymmetry, e.g. 
I51I > I52I, we get 73/72 = {q-l/qf > 1. Therefore the 
three-photon decay channel dominates, resulting in the 
single-photon state generation from any coherent input 
with a large amplitude [19]. For example, for |g| = 10 the 
fidelity of the single-photon state generation can exceed 
99%. For such an asymmetric coupling, in the absence of 
the linear losses, the resulting state is approximated by 
the following one 



[out) 



[g|0)i|l)2-|l)l|0)2] 



(17) 



i.e. the photon is carried mostly by the weakly coupled 
mode 02. The linear losses, however, gradually degrade 
the generated state. For example, for the linear loss rate 
an order of magnitude smaller than 73 , the fidelity of the 
single-photon state generation drops to 86% [19]. The 
effect of the linear loss is smaller for smaller times. How- 
ever, in that case a significant multi-photon contribution 
to the state is still present (Fig. 1(a)). The multi-photon 
component can be a significant obstacle for using the 
generated state in the quantum cryptography protocols, 
where a presence of the two-photon component in the 
signal can be used for security breaking [5]. 

For small linear losses one can obtain a simple estimate 
of the parameters range when there is a significant prob- 
ability of the single-photon state simultaneously with a 
small probability of the multi-photon components in the 
signal. We note that in the Fock state basis only the 
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elements Pn.n+k with a fixed k of the density matrix are 
coupled by Eq. (13). The diagonal elements satisfy the 
following equation 



dpn 



dt 



0'nPn+2,n+2 (18) 



where ct„ = (n + 2)(n + 1)72, /3„ = (n + l)(r + 71273) 
andr„ = n[r+(n—l)72 + (n— 1)^73]. Assuming that the 
mode 62 is initially in the two-photon Fock state, from 

Eq. (18) one easily obtains 

P22 = exp {-2(r + 72 + 73)i}i 

^ 73 + r 

72 + 73 + r/2 
X (exp {-R} - exp {-2(r + 72 + 73)^) (19) 

and Poo = 1 - Pii - p22- 

It is interesting that the simple estimate (19) gives a 
fairly accurate approximation of the single-photon gener- 
ation fidelity (i.e. pn) for times when the initially coher- 
ent state of mode 62 approaches the single-photon level 
(Fig. 1(b)). The approximation is quite precise even 
for the linear decay rate approaching 73. One can see 
that the position of the fidelity maximum is also cap- 
tured quite precisely. Apart from the initial time inter- 
val when the state is far from the single-photon state, the 
estimate (19) describes quite accurately the behavior of 
the two-photon component, P22 (Fig. 1(d)). 

Eq. (19) allows to estimate P22/ Pii at the time when 
the output state attains its maximal fidelity of the single- 
photon component. This latter time, T, is defined by the 
condition dpu/dt = and reads 

1 .l^(l + ll±J^^ 



272 + 273 + F ' ^ ' ^ ' ■ ^^^^ 

This estimate of the fidelity maximum gives a very simple 
result for the ratio 

^ = (21) 

Pii(T) F + 73 ^ ' 

This result seems rather disappointing, since for having a 
small ratio P22/P11 at the maximal single-photon fidelity 
one needs to have the linear loss rate much smaller then 
the nonlinear one, F <C 73. However, if we relax the 
condition of having the maximal possible single-photon 
fidelity, then even a strong linear loss does not eliminate 
the applicability of our scheme for the single-photon state 
generation. Eq. (19) shows that P22 decays faster than 
in the linear case, and pn decays much slower. Even 
for the linear loss rates exceeding 73, our scheme is an 
advantage if compared to the simple linear attenuating 
of a coherent state (Fig. 1(c)). 



The noncleissical states generation with one 
nonlinear and one hnectr side modes 



Now let us consider the asymmetrical case, when only 
one of the side modes is carried by nonlinear (and, con- 
sequently, lossy) media. We assume that K2 = 0, ki = k, 



and the rate of the linear loss is much less than that of 
the nonlinear one, F2 <C Fi. The nonlinear term of the 
Hamiltonian (8), written in the 6-basis, defines the non- 
linear decay channels. Expanding in terms of b\ and 61 
we get the operator coefficients (see Eq. (2)): 

= - (Hw^^^)'^^' = inw^'^^'- 

(22) 

The coupling asymmetry defines the resulting quantum 
state and the output channel (in this case, the local mode 

fli or 02). 

The single-photon state generation.- First, we consider 
the case of gr <^ 1, which leads to the localization 
of the coherent modes as follows bi w 0,2 + qai and 
62 ~ —ai + qa2- Note that, since mode 61 is empty, 
the r.h.s. of Eq. (11) simplifies in this case to the di- 
agonal form with the loss rate F;,^ w Fi -|- q^V2 (note 
that the contribution from the linear losses are negligible 
since q^T2 <C Fi). Resuming, when the modes oq and 
hi are emptied, the leading-order master equation for an 
observable A = A{t, h\,h2) becomes 

^ « i[jy(°'°U]+r6.i?[4]^+722^[(4)']^+73i?[(4)'&2]^, 



with = K{blfbl and the loss rates 



(23) 



72 



AG^/To + Ti 



73 



4GVFo + ri' 



i.e. 72/73 ~ q^ ^ 1, which is the condition for the 
single-photon state generation [18]. This scheme is sim- 
ilar to our previous proposal [19], considered in section 
III A, but involves just one nonlinear side mode. The 
single-photon state is still localized in the nonlinear mode 
(ai). Below we also present the schemes which generate 
quantum states in the linear output modes of the sys- 
tem (whether such a scheme exists for the single-photon 
generation is yet unknown). 

The evolution interval for the generation of single- 
photons is quite similar as in Fig. 1 and given by the 
condition I/73 <C t <C 1/Fi. The sufficient conditions 
(used in the two reductions and the efficiency condition) 
for the single-photon generation read: 



G 

— < 1; K, gFi < — : 
1 1 



G^Fi 



« {qnf, (24) 



where the last condition is needed for efficiency of the 
generator. The consistency of conditions (24) requires 
that FiFo/G^ -C q^, i.e. the efficiency of the generator is 
limited by the losses in the nonlinear Kerr mode. 

The phase state generation. -het us now consider the case 
of (7 3> 1, which also leads to the nonclassical states gen- 
eration. In this case the scheme emulates the two-photon 
absorption, as is seen from Eq. (22). Therefore, for an 
initial coherent state, \z) with |z| ^ 1, the output den- 
sity matrix of mode 62 is very close to the phase state 
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|0) + e^arg(z)|]^^ [15-17]. Note that the output mode is 
still given by 62, however now it is localized in the lin- 
ear mode, b2 ~ a2 — q^^ai. We have the same form of 
the approximate master equation (23) but now with the 



reduced Hamiltonian « {k / q'^){h\Yh2, 73 

and 72 ~ fr§2 7 with 73/72 ~ q^^ <C 1, where we have 
assumed that the linear mode 02 has negligible losses. 
Moreover, since the operating mode 62 is almost local- 
ized in the linear mode 02 the linear loss rate of this 
mode is almost equal to that of mode 02, i.e. we have 
from Eq. (11) in this case: Tb^ ~ r2 + Fi/g^ (the cross- 
term vanishes since the mode 61 is empty). The necessary 
conditions on the scheme parameters in this case read 



G ^ ^ /IT G2 
i V i i 



(25) 



with a similar consistency condition as for the single- 
photon generator, i.e. TiVq/G'^ ^ 1/<Z^- 



was shown that for L = b\ quite different initial states 
decay to the stationary state practically in the same 
time [25]. 

Thus, one can expect that a classical state can be 
transformed into a nonclassical one by the scheme de- 
scribed by Eq. (13) in a short evolution time for which 
the linear losses are still negligible. Indeed, for two non- 
linear side modes, when 73 ^ 72, one should expect a 
significant change of the quantum state to occur for the 
evolution time on the order 



(26) 



Fig. 2(a) confirms this prediction (there the expected 
significant change occurs at 73^ w 0.0004, i.e. exactly at 
the estimate (26)). For such times the average number of 
photons drops by nearly two times. Naturally, the linear 
decay rate should be of several orders of magnitude larger 
than the nonlinear one to play any significant role in this 
case. 



Generation of the multi-photon nonclassical 
states 



Usefulness of the three-mode scheme for generating of 
the nonclassical states is not limited to generation of the 
single-photon and phase states. As it was mentioned in 
the Introduction, the nonlinear decay can be designed to 
produce an arbitrary quantum state. In particular, the 
nonlinear dissipation with the Lindblad generator of the 
type as in Eq. (13), i.e. f{n2)b\, where f{x) is a positive 
definite function, in the absence of the linear losses can 
lead to the stationary superposition of any given number 
of Fock states, or even to the "combing" of the initial 
coherent state by retaining only the components with 
either even or odd photon numbers [24]. 

Moreover, for a sufficiently small evolution time even 
Eq. (13) can lead to strongly nonclassical states. In 
generation of such states our scheme can be much more 
robust with respect to the linear losses than in gener- 
ation of the single-photon state. It is easy to see that 
the modal dynamics in the presence of nonlinear loss can 
be fairly impervious to the linear decay for times much 
smaller than the characteristic time of both linear and 
nonlinear decay. This can be shown by invoking an ar- 
gument from the quantum jumps theory. Let us assume, 
for simplicity, that the evolution is purely dissipative (the 
Hamiltonian is zero). Then the quantum Monte-Carlo 
trajectory, \ip{t)), undergoes random jumps described by 
the Lindblad operator L [21]: 

m))^L\m)- 

The probability of the jump is proportional 
{tp{t)\L^ L\^{t)) . Obviously, the state components 
with the high photon numbers will be "jumping down" 
via the nonlinear decay described by Eq. (13) with much 
larger rate than that of the linear decay. For instance, it 
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FIG. 2. (Color online) (a) Evolution of the photon-number 
distribution for small losses F = O.I73. (b) The photon- 
number distribution at 'ys.t — 0.3; light gray, dark gray and 
black bars correspond to F = 0,0.173 and 73. (c) The corre- 
sponding Mandel parameter Q, given by Eq. (27), for F = 0, 
O.I73, 73 and IOF3 (solid, dashed, dash-dotted and dotted 
curves). We use the solution of Eq. (13) for the initial coher- 
ent state with the amplitude z = 7. Here 72 = O.OI73. 



We also note that even for larger times and a smaller 
average photon number the generated state is quite ro- 
bust with respect to the linear losses. For example, for 
the states with three or four photons on average, the lin- 
ear losses with the rate even comparable with 73 do not 
destroy the nonclassicality (Fig. 2(b,c)). It is interest- 
ing that for moderate linear losses the few-photon states, 
generated this way, can be very interesting. One can ob- 
tain, for instance, a nearly equal mixture of the one and 
two-photon components, Fig. 2(a). 
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The nonclassicality of the generated state can be illus- 
trated by using the Mandcl parameter, 



Q = 



("2> 



- 1, 



(27) 



which is zero for a coherent state. It is seen from Fig. 
2(c) that for the initial coherent state the generated state 
is manifestly a sub-Poissonian one during alH > 0. In the 
absence of linear losses one can distinguish three stages 
of the Q-parameter evolution. First, it sharply drops (for 
the time-scale given by the estimate (26)). Then, when 
the state has a sharp few-photon maximum, the Mandel 
parameter oscillates. Finally, it slowly approaches Qf = 
— 1, i.e. its value for a Fock state. 

For the small and moderate linear loss rate, F < 73, 
only the third stage appears to be affected: the Mandel 
parameter rises slowly towards zero. As long as the linear 
decay rate remains larger then the time-scale given by 
the estimate (26), the first stage is practically unaffected. 
During this stage, the quantum state undergoes strong 
photon-number squeezing. 

Thus the considered scheme can serve as an efficient 
and robust practical source of the sub-Poissonian non- 
classical states. Obviously, all kinds of set-ups mentioned 
above serve for this purpose. For example, a few centime- 
ters of the three-core fiber with the realistic nonlinearity 
and linear losses is enough for reaching Q ^ 0.1 [19]. 

We conclude that the multi-core fiber realizations of 
the scheme described in section II can be quite useful 
for generating of the multi-photon nonclassical states de- 
spite the fact that a strong Kerr nonlinearity is usually 
accompanied by a strong linear absorption. Linear losses 
can be dealt with more efficiently in the scheme with the 
individual damped emitters (quantum dots, atoms, color 
defects etc.) acting as a common loss reservoir for the 
light modes. Moreover, the EIT-likc schemes exhibiting 
a very strong Kerr-nonlinearity [22] might be a promising 
candidate for realization of our setup. 



IV. THE ENTANGLED STATE GENERATION 
BY THE FOUR-MODE COUPLER 

Let us now explore a four-mode generalization of our 
scheme. We assume that a strongly lossy mode Uq is 
coupled to three other modes, one from which is a non- 
linear mode, as, and two other are linear ones, ai,2- We 
consider the Hamiltonian 



Ha = G 



\ geae + ao ^ geal + K{alfal, (28) 



where the coupling coefficients are real and normal- 
ized J2e9e ~ 1- Moreover, for specific localization 
of the coherent modes (see below) the g( are assumed 
to be strongly asymmetric: gi = .92 = £ ^ 1 a-nd 
gs = 1 — 0{e^). The master equation for the observable 



A of the system reads 



-ToV[al]A + r3V[al]A, (29) 



where we neglect the losses in the linear modes ai,2 and 
assume that the strong dissipation of mode ao has the 

shortest time scale, i.e. Fq ^ G, F3,k. To perform the 
steps which lead us to the second reduction, Eq. (7), 
we need to introduce a coherent basis given by the 
unitary matrix U with the constraint Ui^3 = {gi, g2, gs)- 
This matrix can be chosen as follows 








9192 

93 

ff2 



92 

93 



ai 
02 

0.3 



(30) 

In the strongly asymmetric case, the rotation (30) can 
be approximated as bi = — ea-^ + 0(6^), I = 1,2 and 
^3 = 03 + e(ai + 02) + In this case, the adia- 

batic procedure described in Section II eliminates the 
lossy mode (the first stage) and then the coherent 
mode 63 (the second stage). For the evolution times 
t ^ l/(r3 + 4G^/ro), i.e. when mode 63 is already in 
the vacuum state, the two coherent modes 61,2 possess 
two channels of decay, as given by Eq. (7). These chan- 
nels correspond to the two Lindblad generators 
and obtained by the expansion (as in Eq. (2)) of 

Hamiltonian (28) in powers of 63 and 63. In the leading 
order approximation we have: 

if(o.i) = -e^K{b\ + blfibi + 62) + 0{e'), 



hence the corresponding decay rates in Eq. (7) (with 
T = T3+ 4GVFo) read 72 « ^e^Ky{T3 + 4GVFo) 
and 73 « 4e^KV(r3 + 4GVro). The reduced two-mode 
Hamiltonian is iJ(O'O) = e^^K^bl + blf{bi + 62)^- 

The output states of this four-mode scheme are local- 
ized in the linear modes 01,2- Note that the coupling to 
the lossy nonlinear mode 03 induces the coherent linear 
losses with the Lindblad generator b\ + b^ and a reduced 
loss rate F^ = e^Ts, which fact can be easily established 
by expansion of the term FsI'faJ] = r3'D[^j, Uk3b\] sim- 
ilar as in Eq. (11). For below, we can neglect the decay 
channel H^'^'^\ with the decay rate 73, as compared to 
the channel and the linear losses. 

The sufficient conditions for the adiabatic elimination 
of the lossy modes and the condition of the scheme effi- 
ciency read 

4G^ 

ro>G,r3,K; e^K < r3 -h — < CK, (32) 

i 

with the simple unrestrictive consistency condition e ^ 1 
(cf. with that for the three-mode scheme of section III). 
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In Eq. (32) the rightmost inequahty is the requirement 
of the scheme efSciency (the artificial nonlinear loss rate 
72 is stronger than the induced linear loss in 01^2), while 
the middle one together with the left one (already men- 
tioned above) are sufficient conditions for the adiabatic 
elimination of the lossy modes. Under these conditions, 
the strongest nonlinear dissipation, with H^'^'^\ the lin- 
ear losses term TLVlbl+bl] and the reduced Hamiltonian 
together govern the evolution of observables. 
If the liner losses were absent, the nonlinear losses term 
would dominate the Hamiltonian part, since 72 3> e'^K, 
i.e. K ^ + AG^/Tq, which is satisfied due to condi- 
tions (32). To find out what structure the output density 
matrix has, let us simplify the Lindblad generators by 
performing a rotation to the coherent basis: 



c± 



bi ±62 ai ± 02 



V2 



V2 



eV2a3 + Oie^). (33) 



We have iJ^^.o) ^ e^K{3cl - c'i)/2 and H^^-^'^ = 
4e''K(c^)^c^. Then, up to the order the master equa- 
tion (now written for the density matrix, in the coherent 
basis) reads 




FIG. 3. (Color online) The entanglement witness (solid), the 
sum of the negative eigenvalues of the partial transpose of 
the density matrix of modes ai,2 restricted to Fock indices 
n £ {0, ...,5} in each mode, and the corresponding lowest 
negative eigenvalue (dashed). Here the parameters read e = 
0.1, K. = 30, Fo = 10 and F3 = 2 which give 72 = 0.15, 
e^K = 0.012 and Fl = 0.024. We have used the following 
coherent state | — ia, 0.9q:) as the input state with a — 1.6 
(using another input state with the similar average number 
of photons does not change the picture). The figure is an 
average over 17000 quantum trajectories. 



^ = +722?[(34 -^)/2]p + 2ri2?[c+]p. 

(34) 

The nonlinear loss term has the null subspace given by 
the following expansion in the Fock basis of the coherent 
modes 



|null) = ^/(n)5(m)|n,m). 



(35) 



where each n and m take only even or odd values (inde- 
pendcntly), while /(fc) = and 5(0 = ^- Note 

that the coherent state \a, \/3a) belongs to this null sub- 
space. However, only the following Fock basis vectors 



|0,0), |0,1), |1,0), |1,1) 



(36) 



belong to the joint null subspace of and H^^'^\ 

We will further concentrate on the realistic case with 
the conditions (32) satisfied, i.e. when there is also the 
linear loss term and the decay rates are ordered as follows 

73 < e*K< Fl < 72- 

In this case the output density matrix can be found 
only numerically. Since the initial dynamics (i.e. effect- 
ing the two adiabatic reductions) almost leaves the non- 
dissipated modes unaffected, the state of these modes is 
close to the initial state. Hence, one can pass to simu- 
lating the resulting approximate master equation (34) in- 
stead of that for the four coupled modes. While the aver- 
age values can be found numerically without much effort 
for quite large amplitudes of the initial coherent state, the 
entanglement measure (see below) requires knowledge of 
the whole underlying density matrix. This limits the 



efficiency of the numerical simulations and we are able 
to simulate thoroughly only the small amplitude case of 

Fortunately, considering the small amplitude case is 
sufficient for concluding on the general case. This is be- 
cause the nonlinear Lindblad term of the master equation 
(34) admits some integral invariants, which are also in- 
variants of the Hamiltonian dynamics. Indeed, since in 
the absence of the linear losses the photons are removed 
by pairs, the Hilbert space of each of the coherent modes 
is effectively divided into two subspaces, involving only 
either even or odd Fock number states. Hence, intro- 
ducing the identity operators in each parity subspace, 
lev and lodd, the above mentioned integral invariants 
are given by taking the trace of four possible products 
I^^ Ia-2, o'i,2 € {ev;odd}, with the density matrix. For 
t — ^ 00 they define the following diagonal output matrix 
elements poo,oo, P104O: Poi,oi, Piisi- For the coherent 
pure input state |q;i,Q!2) (in the modes c±) we get 



Pnm.nm (^-^) 



where 



dc 



^^lV(")(|ainF(™)(|a2n, (37) 



^^(o)(|an = coshdap), F^^\\af) = sinhdap). 

We have verified numerically that for jap > 2 the two 
functions F'-^^x) differ by less than 4%, i.e. all Pnm,nm 
are very close to each other for such arguments. Noticing 
that the output density matrix is close to a pure state 
(and not proportional to the unit matrix) we conclude 
that the output state is an entangled state for any initial 
coherent state and small linear losses. 
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V = (a) 


Jllhlllllll 


llllllllljlll...... 



001020304050 01 1121 3141 51 0212223242520313233343 53 0414243444540515 25354555 




001020304050 01 11 21 3141 51 0212223242520313233343 53 0414243444540515 25354555 

Fock basis indices; n,,n^ 



FIG. 4. (Color online) We give the structure of the eigenstate 
of the density matrix corresponding to the maximal eigenvalue 
(i.e. the maximal probability) of the density matrix for three 
different propagation times, 72t = (a), ~ 5.95 (b) and 
72t = 18 (c). Here the abscissa is the ordered Fock basis 
indices nin2 of the state amplitude. The parameters are as 
in Fig. 3. 



The predictions were checked with the quantum 
Monte-Carlo simulations of the master equation (34) in 
the local basis ai,2 using the "quantum jumps" method 
[21]. The numerically obtained output density matrix 
indeed describes an entangled state of the two linear lo- 
cal modes ai and 02, see Fig. 3. We use the negativ- 
ity as a measure of the entanglement [26, 27], but with 
the minus sign in front (for convenience of comparison 
with the smallest negative eigenvalue) . It is equal to the 
sum of all negative eigenvalues of the partially transposed 
with respect to the local modes ai^2 density matrix, p^^: 
■^{p) = J2e<o ^(P^^) i'^^^ *° technical reasons, in cal- 
culation of the negativity the expansion over the Fock 
states was truncated at n = 5, however, this turns out to 
be sufficient to capture negative eigenvalues [28], more- 
over, the higher Fock number states are populated sig- 
nificantly only for the relatively short time). The output 
density matrix turns out to be close to a pure state (for 
a pure input state): the largest eigenvalue of the output 



density matrix drops during the evolution only to 0.85 at 
72 1 — 18 in Fig. 3. The significant elements of the cor- 
responding eigenstate of the density matrix involve the 
trivial null states given by Eq. (36) and few Fock states 
with higher population of the local modes ai_2, see Fig. 
4. Our simulations show that this seems to be the only 
outcome for the output density matrix at least for small 
amplitudes of the coherent initial states. In view of the 
above arguments, this result must be expected also for 
arbitrary amplitudes. 

V. CONCLUSION 

We have shown that it is possible to generate the non- 
classical states of light by emulation of the nonlinear ab- 
sorptions using the usual Kerr-nonlinearity, linear mode 
coupling and the strong linear losses. We have concen- 
trated on the potential capabilities of the scheme to gen- 
erate various quantum states in the simplest possible se- 
tups of three and four coupled modes. We have found 
that, additionally to the single-photon states, also the 
phase states and the entangled states can be generated 
by our method. It is interesting to observe, that though 
our scheme depends significantly on the Kerr nonlinear- 
ity (which must be present at least in one of the modes) 
the output states can be localized in the linear modes, as 
in the case of the phase states and the two-mode entan- 
gled states. Though, we have not been able to present a 
similar scheme for the single photon states, there is some 
hope that such a scheme may also exist. Having the 
output state localized in the linear modes of the system 
prevents its subsequent degradation by the destructive 
linear losses, accompanying any known strong Kerr non- 
linearity, but quite small in the linear media. We believe 
that our scheme can be realized in a practical setup with 
the current technology. Some perspective candidates are 
the optical fiber coupler or the multi-core fiber, the elec- 
tromagnetically induced transparency media, where the 
"giant" Kerr nonlinearity has been demonstrated [22] . 
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